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ABSTRACT 


THE DYNAMIC PROPERTIES OF THE FULL JOURNAL BEARING, 
CONSIDERING CAVITATION 


by Richard A. Kaufer 


Submitted to the Department of Naval Architecture 

and Marine Engineering and the Department of Mech- 

anical Engineering on 20 May, 1966 in partial 

fulfillment of the requirements for the degree of 

Naval Engineer and the degree of Master of Science 
in Mechanical Engineering. 


The unsteady motion of a circumferentially grooved journal 
bearing is analyzed and programmed on a digital computer, using 
the short bearing approximation of Reynolds! equation. The 
boundary conditions of continuity are used, and particular em- 
phasis is placed on determining the effect of the size and 


shape of the cavitation zone on the unsteady characteristics of 
the bearing. 


Tabulated data consisting of bearing loads, load angles, 
and spring''and damping coefficients are presented for a wide 
range of bearing dimensions and operating conditions. A com- 
puter program is also written which can be used to obtain a 
wider range of values than are tabulated, and/or as a sub- 
routine in an analysis of a journal bearing supported systen,. 


The analysis shows that the dynamic properties of a jour- 
nal bearing are greatly influenced by the shape and extent of 
the cavitation zone, thus indicating a need for experimental 
work with the primary object of ascertaining film behavior 
under unsteady conditions. 


Thesis Supervisor: 


Herbert H. Richardson 
Associate Professor of Mechanical Engineering 
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NOMENCLATURE 


damping coefficient, as in By,; lb-sec/in. 
Bman C 


dimensionless damping coefficient, 
Pa LR 


radial bearing clearance; in. 
bearing center 

journal center 

perturbed journal center 
perturbation of journal center; in. 
journal diameter; in. 

journal eccentricity; in. 

total bearing load; lb. 


total steady state bearing load; lb. 
25, 
Pa LR 
bearing forces in 1 and 2 directions respectively; 1b. 


dimensionless bearing load, 


Steady state bearing forces in the 1 and 2 directions 
respectively; lb. 


bearing forces in r and s directions respectively; lb. 


radial bearing clearance at any point; c(1l+ncos 60); 
in. 


radial bearing clearance along the upstream boundary 
of the cavitated zone; c(l+ncos 0*); in. 


fovaz 


spring coefficient, as in K,,; lb/in. 
K pane 
Pp, LR 


total bearing length, excluding circumferential 
groove; in. 


dimensionless spring coefficient, 


eccentricity ratio, - 
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pressure at any point; psia. 
feed pressure; psia. 
feed pressure, (Pp, -P,) psig. 
ambient pressure; psia,. 
Po 
pressure ratio, — 
Pq 
yolumetric flow in the m direction; in*/sec. 
journal radius 


radial and tangential reference axis 


fluid velocities in the x,y, and Zz directions res- 
pectively 


perturbation of journal center; in/sec. 


film width at any point measured from the circumfer- 
ential groove; in. 


horizontal and vertical reference axis 


film coordinate axis 


a 
dt 


width of a oil striation at the downstream boundary 
of the cavitated zone; in. 


width of a oil striation at the upstream boundary of 
the cavitated gone; in. 


angle between bearing load and line of centers; ra- 
dians 


rotational speed of journal; radians/sec. 


relative tangential velocity between journal and 
bearing, (W- 2%); radians/sec. 


dimensionless squeeze, AA (Lya 
Pa c 
viscosity; lb-sec/in*. 


angle between line of centers and any point on the 
journal; radians 


shear stress; psia. 
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INTRODUCTION 


Background 
The knowledge of the dynamic properties of a journal bear- 


ing is necessary in the design stage both to ascertain the 
inherent stability of the bearing itself and to predict the 
overall response of a bearing supported system to external and 
4nternal excitations. Some of the many problems which can be 
encountered are self excited half speed whirl, excessive am- 
plitude of vibration of the journal within the bearing near 
the critical speed of the system, and the inability of the 
system to recover from a dynamic load. 

Since Reynolds‘derived the three dimensional, continuity, 
differential equation for determining the pressure distribu- 
tion in an oil film between two moving surfaces, all journal 
bearing investigations have used Reynolds! equation, usually 
in one of two simplified forms because of its mathematical 
formidability. In the two dimensional or infinite bearing sim- 
plification, the pressure induced flow in the axial direction 
is neglected. In the short bearing simplification, the pres- 
sure induced flow in the circumferential direction is neglected 
by assuming that it is much less than the velocity induced flow 
caused by the relative motion between the journal and the bear- 
ing. | 

However, defining and applying the proper boundary con- 
ditions has ‘proved to be the most difficult task in obtaining 
& good solution to Reynolds! equation, even in one of the re- 
duced forms and under steady state conditions. Sommerfeld, in 
1904, was the first to obtain an exact solution to Reynolds! 
two dimensional equation, but his solution allowed negative 
pressures in the diverging portion of the bearing. This led 
to the theory that under certain conditions of loading, vapor- 
ization could occur because a fluid’ is unable to support any 
appreciable tensile force. This theory has since been verified 
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by making observations of tranlucent bearings which clearly 
show film rupture[1],[2],[3],[4]. Theoretically, the cavita- 
tion pressure should be the vapor pressure of the fluid, but 
4n actual practice the cavitation pressure has been found to 
be very nearly ambient pressure because of entrapped air and 
other gases. 

In 1952, Ockvirk[5] developed a steady state solution of 
Reynolds! equation using the short bearing approximation and 
assuming that the entire unloaded half of the bearing was ca- 
vitated and at zero pressure. More recently, Fedor[6] and 
Donaldson{[ 7] obtained solutions of Reynolds! three dimensional 
equation using series solutions and applying Sommerfeld's 
boundary conditions, but requiring that the feed pressure be 
sufficiently high so as to suppress cavitation. However, this 
is an unrealistic situation in most bearing applications, par- 
ticularly those in which stability is a problen. 

Two recent analyses of unsteady motion are of special 
interest in this investigation. Holmes[8] extended Ockvirk's 
solution to the particular case of a rigid shaft supported by 
Symmetrical bearings. Richardson[9] also used the short bear- 
ing approximation, but defined the cavitation boundary more 
rigorously by assuming that there was a zero pressure gradient 
along the boundary and by assuming that the shape and extent of 
the cavity was a function of feed pressure and squeeze..Rich- 
ardson obtains a solution to his analysis by digital methods. 
However, both Holmes and Richardson violate continuity since 
they do not consider the oil entering and subsequently leaving 
the cavitation zone. | 

In 1957, theories which took into account continuity of 
flow through the cavitation zone were independently developed 
by Floberg[10] and Rightmire[11]. Experimental work by Flo- 
berg and Ishii[4] has indicated a close agreement with the 
continuity theory under steady state conditions. 
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Purpose 
The purpose of this investigation is to analyze the un- 


steady motion of a circular, circumferentially grooved journal 
bearing, using the short bearing approximation and the bound- 
ary conditions of continuity. A circumferentially grooved 
journal bearing was chosen because its line source makes the 
feeding arrangement independent of the journal attitude within 
the bearing. The results of this analysis will be presented 
in two forms: a computer program which can be used as a sub-= 
routine in a systems analysis, and in the form of tabulated 
data.of spring and damping coeffecients for use in hand calcu- 
lations. 


ANALYSIS OF UNSTEADY MOTION 


Assumptions 
In arriving at a solution of Reynolds! equation which 
will adequately describe the unsteady motion and dynamic pro- 
perties of a circumferentially grooved journal bearing, the 
following assumptions will be made: 
1. The fluid is Newtonian and incompressible. 
ee tne viscosity jis constant. 
3. The flow is laminar. 
4, Tne mass of the fluid is negligible. 
5. Surface tension is negligible. 
6. The fluid adheres to the journal and bearing surfaces. 
7. Lnere is°ne distortion of the journal or bearing 
Surfaces due to hydrodynamae-pressure: only ignt 
and moderately loaded bearings will be considered, 
See lne TLvid.Titim 1S socthin thav it has no curvature. 
Seeithe fluid film.is so thin that there is-~no pressure 
gradient across the film. 

10. The vapor pressure of the fluid is equal to ambient 
pressure. 

11, The pressure induced flow in the circumferential di- 
rection is negligible compared to the velocity 
induced’ flow in the circumferential direction 

12. In the cavitated region, the fluid divides into an 
infinite number of streamers. 

The last two assumptions will be discussed in greater 
detail in later sections; the remaining assumptions are self 
explanatory. 


Journal Bearing Oil Film Theory 

Figure I illustrates the general case of the unsteady 
motion of a’ journal within a fixed rigid bearing. At any 
instant in time the journal has a rotational speed W about 
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Figure I 


Journal-bearing motion referred to a fixed point in space 
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Figure If 


Journal-bearing motion referred to the line of centers 


Figure III 
Oil film forces in the x direction 
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Figure IV 
Oil film forces in the z direction 
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Mos. Own center CJ, an orbital speed of g about the bearing 
center CB and a radial velocity e€. However in analyzing this 
Syocem, it is more convenient to fix the line of centers: in 
Space and use that line as a reference as shown in Figure II. 
Then the journal has a clockwise rotational speed of (W - B) 
and the bearing has counterclockwise rotational speed of D. 
The forces acting in the x and z directions on an ele~ 
mental cube dx dy dz are shown in Figures III and IV, where 
“ee ond w. are the fluid velocities. inthe x, y, and Z di- 
Tections respectively at any point. Equilibrium of forces 
mn tne x and y direction yields: 


OP _ Ox, Cae 
gx dy 32 (1) 


pire ots, Ola (2) 


oe By ax 


With the assumption that the film is so thin that has no 
pressure gradient across it; 


oP _ 
ay 


ie) 


The sheer stresses are obtained from Newton's law for laminar 
viscous flow[12]: 


Substituting these expressions into equations 1 and 2 gives: 
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Again considering the fact that the film is very thin in the 
weoebection, so that the velocity gradients in the y direc. 
tion are much greater than the velocity gradients in the x 


and Z directions: 


Bp 1, 0° u 
pc ye (5) 
oF Hy 8 (6) 


Withthe assumption of constant viscosity, equations 5 and 6 
can be integrated twice: 


aD 
3P 
Gato soe¥ 71 taCs¥ut Ca (8) 
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where C,, C,, Cs, and C4 are constants of integration. 
From Figure II, the boundary conditions are: 
at y (W- Z)R, w= 0 


at y=h: u= - QR, w= 0 
where h is the film thickness at any x. 


I 
(o) 
S 
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Applying these boundary conditions to equations 7 and 8: 


eae) Ty! aD 2_ _ WR \ ty & 
6 nla re 


With the assumption that the pressure induced flow in the 
circumferential direction is much less than the velocity in- 
duced flow due to the relative motion between the journal 
and the bearing, equation 9 becomes: 


i ~ “Say Ha(US Eg) R (fay) 


The boundary velocities in the y direction are: 


Se 


cee Ope ae = 


at y eae ieee 8) 


By acontinui ty: 


Substituting equations 10 and 11 into the continuity equation: 
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Integrating the above expression from y = 0 to y = h and 
applying the boundary velocities in the y direction, the 
final result is the short bearing approximation of Reynolds! 
equation for unsteady motion: 


Supe | 64 Mae Sete 
gz? ne bx ne fe 


where A. = (W- 29) 


This equation is much easier to manipulate and solve 
than Reynolds! complete equation, and for moderately loaded 
bearings with length to diameter. ratios of less than % the 
loss of accuracy introduced by the short bearing simplifi- 
cation can be expected \to be of the same order of magnitude 
as some of the other assumptions which were necessary to 
make. A comparison between an exact solution obtained by 
Donaldson [7] under steady state conditions (= 0, é = 0) 
and the results of this investigation is presented in 
Appendix A. It is proposed that the same deviations could 
be expected in the case of unsteady motion between exact 
and approximate solutions. 

Again assuming constant viscosity, equation 12 can be 
integrated twice with respect to z to obtain the pressure 
distribution on the bearing surface: 


y anO7t 
na = — 75y 2 + C24 Ce (13) 


10 


where C, and Cz are constants of integration to be de- 


termined by the film boundary conditions. 


Expressions for the unit flows in the x and z directions 
will be required in the next section for determining the cav- 
tation boundaries. These expressions can readily be obtained 
by integrating the velocities u and v across the film as 
follows: 


h 
a, ={ vay 
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Bupstituting in”equations LO and 11 the final result_is: 
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Film Boundary Conditions 


The unwrapped oil film for one half of a circumferen- 
tially grooved journal bearing is shown in Figure V. The 
approximate shape of the cavitation zone as observed experi- 
mentally [1], [2], [3], [4], 1s defined by the boundary 
1 - 2 - 3; for steady motion, point 1 occurs immediately 
downstream from the point of minimum film thickness (x =nR) 
where the rate of film expansion is the greatest, but in the 
case of unsteady motion it must be assumed that cavitation 
could commence at any x. Since the vapor pressure of the 
fluid is assumed to be ambient, then the entire cavity must 
be at ambient pressure, with the result that one of the 


it 
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boundary conditions is: 


where w is the film width at any x and is equal to L 
where cavitation is not present. 


Applying this boundary condition to equation 13, the result is: 


(2% Swe) +6 pz, = Ww) (16) 


PLeune. 
The unwrapped oil film 


| | \ cavitated 
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With the assumption that the fluid adheres to the sur- 
faces of the journal and bearing, the fluid must split into 
one or more striations upon reaching the boundary 1 - 2 in 
order to satisfy continuity and to prevent subambient pres- 
sures from occurring. This is shown in Figure VI for the 
particular case where h,>h, . 


a2 


Figure VI 
Continuity of flow across the cavitation boundary 
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The upstream boundary 1 - 2 is determined by considering 
‘the flow through an elemental volume h dx dz at A. Then the 
‘flow entering the boundary at A is 


QO =O dx + O.dz 


bN 
From equations 14 and 15 this becomes: 


. hn? ap 
= 4 = - —-— d 
Roe s(W- 28)Rhdx ise Zz 
Since there are no pressure gradients in the cavitated Zone, 
the only flow is velocity induced flow caused by the rela- 
tive motion between journal and bearing. Thus the flow 


leaving the boundary at A must be: 


Q ouT 


= $(W - 2%)Rh dx 


Therefore, by continuity of flow across the boundary: 
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or: 


along boundary 1 - 2. 


Applying this boundary condition to equation 16 to evaluate 
the constant C, the pressure distribution in the region 
bounded by z = 0 and 3 = 4 -1 = 2 is found to be: 

SE 


h 
p-p Lindh 


Geta s ee 2y)(z*-2zw) + p'(1 - =) (@atiz3)) 


where p!' = ir 1 


The first term on the right side of equation 17 is the 
hydrodynamic pressure distribution and the second term is the 
hydrostatic pressure distribution. The slopes of these two 
terms at Z = WwW, are, respectively: 


‘ h 2 
ta gs pyc de gs me a0). Cie) 
B= hid lal ax 
= pt 
pe ae ee (19) 
| a hs a 


The existence of cavitation at any x upstream from point 2 and 
downstream from point 3 can be ascertained by letting w = L in 


the above expressions. Then the condition for cavitation to 
Gecur is: 


3LA 
h? 


Sito ep eg 
= =2y) 7 (=) (20) 


(RA 
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i.e. the positive slope of the hydrodynamic pressure is 
greater than the negative slope of the hydrostatic pressure. 
This is illustrated in the left diagram of Figure VII. 


Figure VII 
Axial pressure distribution in region l - 2 


hydrostatic 


i 


equation 20 (unstable ) equation 22 (stable) 


If the condition for cavitation is satisfied(equation 20) then 
w adjusts so tnat: 


BL (H-° 


Substituting equations 18 and 19 into the above expression, 
the film width along boundary 1 - 2 is 
3 
2 pth 


des ic (21) 
oh 5° | 
3M (RASS ~27) 
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Making the substitutions: 


h cn 
ig 77 sine 
y =- e cos 9 


3A/ (cnAsin®©O- 26 cos 6) (225) 


If equation 20 is not satisfied, then: 
Ce aa OF 


It can be seen that equation 22 is invalid for the down- 
stream boundary 2 - 3 by considering the steady state case, 
Then equation 22 becomes: 


-p'h?® 
W= [3747 cnW sin © 


Under these conditions, the expression under the square 
root becomes negative at angles greater than 9 = 2n, implying 
that the cavitation zone cannot extend past 9 = 2m. Physically 
this would not necessarily be true for high journal speeds and 
low feed pressures. 

The downstream boundary is determined by considering con- 
tinuity of flow through the cavity as illustrated in Figure. 
VIII. Let the width of one striation at the upstream bound- 
ary be Az. Then the amount of fluid in the striation as was 
previously determined is: 


$3ARh* Az 
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where h* is the film thickness at Zz = w on the upstream 
boundary. 


Figure VIII 
Continuity of flow across the downstream boundary 


1 


With the assumption that the fluid adheres to the journ- 
al and bearing surfaces, then by continuity the same amount of 
fluid must leave the cavity at Z = w at point B on the down- 
stream boundary. Since we have assumed that there are no 


pressure gradients in the cavity, there can be no motion of 
the striation in the z direction. 
Summing up the total flow in the unit volume h Ax Az 


at B, the result is : 


a { 
ZARH*A z - a 32 | Ax = /ARh Az (24) 
ly 
From equation 17: 
oP | 344 gn ® jy 
22 = (RAS - 2y)w - = (25) 
W ‘ ‘ 


evan - 


sae feet &: 
oe oie oes G2 Wal; 


a 


* 


Substituting equation 25 into 24: 


Ane ae 22 2 (RAS2 . yyw! 
2 ce ee ons (R Neue og woos en ae 


Assuming an infinite number of striations, and then taking the 
limit as Qx-edx and Az—daz: 


a] 
ZARh* dz -| Faas - 27) - 22 jex- 2 ARh dz (26) 


Making the substitutions: 


dx = Rdée 
ae = dw 
|W 
ca = - S sine 


y= -e cose 


Equation 26 becomes: 


if c . hép'| 
3ARh*dw -| dw(- cnAsin@+ 2€ cos 6) - ——iR de 
Le Wi 
= $/\Bh dw 
Ox s h?p! 
dw gw(onAsin 6 - 26 cos 6) + 64w Bek 
aa A(h = h*) 


The film width, w, in equation 27 is in terms of two vari- 
ables, dw/d® and h*, and so must be solved by an iteration 
method. 
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CF) 


O11 Film Force 


.The oil film forces which will be considred are: 
1. Hydrodynamic, hydrostatic and shear in the complete 
film region bounded by. z = 0 and 4-1-+2-+3 +4 
Mm Fi gure: Vs 
2. Shear in the cavitation zone. 
The differential hydrodynamic and hydrostatic force on a 
strip of the journal w, Rdé@is: 


oe 
ak, =F) paz kas 


Z=0 
mupstituting in p from equation 17: 
ue h " z 
dF, [| Mat ~ 2y)(z? - zw) + p'(1 - 2)lazRae 
bay hs x W | 
Bo ae ~ 


which when integrated becomes: 


Resolving this differential force in the 1 and 2 directions as 
shown in Figure IX: . 


dF, = dF, cos @ 


aR 


P2 = dF, sin 6G 


The tangential shear stress in the complete film region 
due to the relative sliding velocity is (see Figure II): 


= we - Gecos® + ésiné) 
However, #e<<WR and éK«WR, so that: 


EY) 
Ga = 


9 


4 


oS & i + AY 
‘ t 


Beets S-Otem 2 ete a? eer? ete +: b crt garek 
fount A0Xe 


oe the a 


Figure IX 
Coordinate system for oil film forces 


Then the differential force due to shear on a strip w, Rdé@is: 


Substituting in the value for U: 


MUN Rea 
dF... - Lee 
= nh 


Resolving the shear force in the 1 and 2 directions: 
aFy, = - dRsinoede 


dF, = dF cos 6 dae 


7 
In the cavitated zone, there is no hydrostatic pressure, 
Since p = p 


a 
there is no pressure flow and thus no hydrodynamic pressure 


- Also, since there are no pressure gradients, 


Smt s cS . 
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The shear forces in the cavitated zone are obtained by 
assuming an infinite number of striations as in the preceding 
section. Refer to Figure X where one striation is shown. 


Figure X 
Continuity of a single oil striation 
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By continuity: 
SARN*S* = ARNG 
where &* and 4 are the striation widths at the upstream 


and downstream boundaries respectively. 
Then: 


or: 


The differential shear force due to one striation at any @ is: 


df, =UsR de 


ak 


Substituting in the expression for & from above: 


AWR h*o* 


dt; n = R de 


In the limit, as S5*—>dz: 


ry 2 
df, = ate h*dz ae 


Then the differential force due to shear across the complete 
width of the cavitated zone at any 9 is: 


L 
Ad 2 
dF. Lf: WH" n* az de 
7c} h 
£z=w 


where the integral must be evaluated numerically. 


Resolving the cavitated zone shear force in the 1 and 2 direc- 
tions: 


aFe =-d. sine de 
aF. = dF. cos iS) ae 
“c2 7c 


Summing up the ferces in the 1 and 2 directions and making the 
substitutions: 


y = - é cos 6 


Then the differential forces in the 1 and 2 directions are: 
{wR 


anaes WRP 
| one 


q 
cos @ 


(J.cnsins - 2é cos @)cos 6+ 


A 2 bn4 } 2 
_ MUR Bean _ MWR ifs 


sam. |'d6 55 (25) 
h ne | 
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—— (Acnsin® - 2€ cos 6@)sine+ 


E 
where I -| h*daz 
WwW 


Sr iT 
=/ dF, 


7 
F =| ak 
ol 2 
Jo 


which must be evaluated numerically. 


MWR? w M( (URI e 
oi eGR aes © + See oe (29) 


Then: 


TZ 
Fr 
‘ 


The attitude angle is 


be 1,2 
g = tan Se (30) 


Spring and Damping Coeffecients 


Refer to Figure XI. At a given instant of time the jour- 
nal center is at CJ and is being acted upon by the forces oe 
and F, which were evaluated in the previous section. Now the 
journal is given a perturbed displacement, d, and a perturbed 
velocity, v, to the position CJ'; then the forces acting upon 
the journal are of the form: 


Pt Re (Sh) ety BiG, S.) 


1d 


He 


° 


Hop HKe(r,s) d+ Bo (r,s) 


where F,, and Fz, are the steady state forces at CJ; and K 


and B are spring and damping coefficients. 
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Figure XI 
Perturbation of journal center 
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In evaluating the previous equations it is first nec- 
€ssary to. transfer the forces F, and Fz to the r and s axis. 
However, it is more convenient to continue using differential 
forces; thus: 


dF, cosx~ dF, sina 


dF. dF,sinA+ dF,cosa 
Since a small perturbation is being considered: 


Coc". << CB CJ! 


Then : 


24 


; | a : a | ; 7 4 
oe Z : > * ) allt 
Elen ». ite ivr og -% Gad 
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7 is mote Nh «one 1 ace 
. > 
= .> 4 


Sine, "t + dug te 4d - 2a 


these jen gale? rt tniktac ie t8q tioes | 


1S. ao op ee 


ai 
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With these approximations, the differential radial and tan- 


gential forces become: 


Substituting the eauations 28 and 29 into equations 31 and 32: 


wR , 
dF, = ae cos 6 | - 2e cos 6 +/\nc sine | +2 pwRoos 6 


df 2 AA LIN — 
Ss Was MW RFT (Hat aE 


no - —=—sinoe 
ig ie 
fe 7 
ba 2ecose6 ++ nesin6| + $p'wR sine 


\ 


——— cos 3) de 


ene i 


c 


( s (MWR iar IN op 
dF, = ) on cos 6} -2ecos6 + ensin6/|+ 2 p'wrcos 6 


{/1) D2 2T 
elie taht rea ey enti sag) erg 
a n® 2n° 


peer 


- 2ecose9 + +incsinG! + = p'wR sin ® 


i 
ES 


| HW RS: UUWR? F ) 
+ detneatelle 4 6 + ee a 
h h? 
Let: 
é=r 
h = c(1 + ncos @) 
: S 
G = on 


Substituting the above into the equations for ak. and dF. 
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and rearranging: 


. CY wR 2rcos 6 WWne sine 
Gh = <———cos 6 | = Syl Oso ee ANG 
r 205 il meos'6)* «(itn cos. 6) 
Past Oe ee ee Bow sin © 
(1 +ncos Bye ats (l+ncos 6)c 
Va ae L sin 6 Mw®Rsine fF 2rs 
c*(1l+ncos 06)? 20° | en(l1+ncos 6)* 
Us. sine 2ss sin® } pwRs sin 6 
(l1+ncos 6)? ~ en(l+ncos@)* | 2cn 
AW R®wscos6 AWR*Iscosé6 ) 
~ Sa de 033) 
e*n(l+ncos6 c®%n(l+ncos e6F 
tee M w*R cos 8 2rs cos @ ius sine 
cae Ze (l+ncos 6)%cn (1 +ncos @)* 
2ss sin® i Pswhcos6 WHiW)R*ws sin 6 
en(l1+ncos 6)* | 2cn e7n(1l+ncos 6) 
A{W R®Is sin © Adw®R sine f. 2r cos © 
e°n(1+ncos 6)? 2c% | (1+ncos 6)® 
UW) ne sin é 2s sin © : 
—________— - ——________ | + 3 ywisine 
(l+ncos 6)° (megs 6)? 5) 
MW Rew eos 6 MWR®tcos 6 
ey ee ae (34) 
c(1+ncos © c?(1+ncos 6)? 


Equations 33 and 34 are non linear in n. However, with the 


assumption that the perturbation is small, the equations can 
be linearized. From Figure XI: 


en, +r = cncosnv sen 


OZ: 
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where = << n,and is time dependent. 
By Taylor's theorem: 

f (n+ =) = f(n,) + = f(n,) 

neglecting second order and higher powers of = . 
Then, by also neglecting other second order terms: 


s Eats S 
oy E(no +=) WY > f (n,) 


5 f(n, +=) =~ rf(n,) 


r Toe Tr 

eto tee ne 5 f(n, ) 

s f(n,+~) ~% s f(n,) 
ie) fe Ueteeed 0 

rs= 0) 

ss=0 


Linearizing the terms of equations 33 and 34 as above: 


Ad wR 2r cos e iy onsin e 
20% Ia Cle+ neos.0)F a (d+ m.cos' 6) 


Wrsin® (1+ ncos®© - 3ncosée6 25 sine 
2 See = = SSS 
1+ncos 6)4 (1 +n cos 8)? 


ANWR? w sin © is AWW R7wr sin®Ocos © 


+ gp'WRecos® - 


e(l+ncos @) c*(1+ncos 0)? 
. MWR*T sine 2MWR*Irn sin 6 cos 6 
c?(1+ncos 0)? c*(1+ncose)® 
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= 


i aprade 64 eC ita vag afctiasps Ww cine? syd 


ne ti) ce y a 
7 sgt est -] si 


M w*RsWsin? 6 PwRssinOo AMWR®wscos 6 
2c°7(1+ncos 6)? 2cn c?7n 1 +ncos 6) 


Ad) RF s 6 
" = Is co : ae (35) 
e°n 1+ncos 6) 
ara a Alw?RsWoos 6 sin © PswRcos@® Ai(yvR ws sin® 
S  ) 2e%(1+ncos 6)3 2cen c®=n (1+ncos 0) 
MW R?Is sin® JA4 w*R sin © [ 2r cos © 
c*n (l+ncos 6? 203 {| (1+ncos 6)* 
iS Wne sine Wer sin @(l+ncos6 - 3ncos 6) 
(1+ncos 6) (l1+ncose)‘ 
3sine 7 Li) B2wcos © 
eee OS a : Perera etme 
(l-+ncos 6)" | c(l+ncos 6) 
MUWR*wr cos? 6 MUR?I cos © 
~ 92(1+ncos ©)? c?(1+ncos 6)? 
2MWR®Ir cos? 6) 
- FEST fT TE Oe (36) 
c*(1+ncos 6) \ 
Equations 35 and 36 are of the form: 
dF, = - dK,,r - 4B,,T - dK,,8 + 43,8 + dF, 
7, = = Gi dB. c e. 1B..# + ak 
GF, 2 = GK,.8 = GBeg 8 + GKp,F + GBs. F + Pay 
where: 
Af (Uw? R 
Sin COs oll +n COs 6 = 3n cos 6) 
dK» j= 5 a 
MWR*w sin © cos 6 2AWR*I sin @cos 6 a6 (37) 
c#(1+ncos 6)* c?(1+ncos 6)* 
3 2 
pee ee oe ae (38) 
e“{1+:ncos 6)* 


28 


; | 
«=&nrpec + Ds ents =—— aA - 
a (a eoen + fT" 


— 


Sy ewe BL sin” 0 PiwR sin® | AW R3wcos @ 
rs) 2c? (1 + ncos e)° 2cn c7n (l+ncos @) 


MW R®I cos @ 
pee 108.831 f (39) 
pe oh 8 Sa a ato KSA 
Ad w°R cos 6 sin 6 
Sh SS AO ee Me 
rs e*(1+ncos @)* (40) 
ane Mw RW cos@sin® pwkRcos ce AAW R?7w sin 6 
Ss 7 2c7(1+ncos e)* 2cen c7n (1 +ncos @) 
)R°I sine 
Ahi ee (41) 


c°n (1+ ncos 6)? 


M4 wR sin? © 
dB... = Soe (42 
ss c°(1+ncos @6)* ) 


AV w*®R sin? © 


ak... = - Ll+ncos9 = 3ncos @ 
ee 2c%*(1+ncos @) ( S 
J 2 2 2ZHMWR*I cos? © 
eee wcos” © Ueki UJ os Bae (43) 
o7(1+#ncos 6)" c®(1+ncos @) 


a Al w°R sin @cos © Se (nin) 
sr Fs 6 (1 +ncos 6)? 


then: 
oat 
Ade = | AK jy re iL) 
oO 
ai 
se = f GB nnn 
o 


which must be evaluated numerically. 
The final result is: 


ees Me Ps sree es 
Py Kt Bt K.. 8 + B.S + F, (45) 
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Fos -'K.s Pee oe dete oe te be Oe Fe (46) 


However, it is generally more convenient to work with 
cartesian rather than polar coordinates. Therefore the in- 
tegrated forces, and spring and damping coefficients will be 
transferred to the x - y axis as shown in Figure XII. 


Figure XII 
Journal motion in x = y coordinates 


CB 


Sy. 
From Figure XII: Vx 
r= ysing +xcosg (47) 
s = ycosg -xsing (48) 
r= ysing +x cos B (49) 
s = ycosG -xsing | (50) 
F, = F cos §- F; sing (51) 
R= F,sin %+F;cos g (52) 
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Equations 47 through 52 are substituted into equations 45 and 
46; 


hy 
U 
' 


Kp, (¥ sing + xcos Q@) « By-(y sing + xcos )+ Fi, 


Kps (y cos £ -xsinf)+ Bp.(ycos# -xsing) (53) 


Fe @ = Keo (yoos@ - x sing) - Be, (y cos @ - xsinf)+F, 


“ 


+ 


K.,»(y sing + xcos$)+ Bsy (y sing + xcos@ ) (54) 
Substituting equations 53 and 54 into equations 51 and 52: 
cat app sing cos fy - Ky,cos® f x - By, sing cos Py 
- By, cos? J x - K,,cos* fy + Ky, sin cos $x 
+ B;, cos? g ne By, sing cos Px + FF, cos. 9 
+ K,, cos@ sing y - Kgs sin? 9 x +B,s,cosf sing y 


- Bes sin? g x - K., sin® g y - Ks, sin @ cos gx 


- B,, sin? g y - B,, sing cos#x - F 


sing 


ar 
“Ss 


2° 
F, = - K, sin? g y - K,,, sing cos $x - B,, sin? gf y 
- B,, sing cos x - K., sin cos fy + Kiss sin® @ x 
+B sing cosfy -B,_ sin? fx + F,, sin @ 
- K,, cos*  y + K,, sing cos fx - B,, cos*f y 
+ By, sing cos Gx + K,, sing cos Py + Keg cos* g x 


tbo. sin cos fy tobe. COS. x + pas cos g 


31. 


fa see's, 4 


a 7 i 
: Par ® 05/8 wii, = he ogee 
> 7 


aase ag! = BE toe Wh a56. - Ea. rn 


; Pinte ate =e oe Oats wh = % . "alt ‘ 
anny 

om, = it @Aises Male g Bin in me Gore oe 

= : . 

Reins nvenya tame Gots at #. 
iP, lye - - Amr 


ee8 iat 
itt co a rt ‘ 


Then: 


PKK By is Kv Seyy com re G55) 

DP age SS phil a ab Ml a (56) 
Where: 

K,, + Kp, cos*~ + (K,,- Ky),) sing cos f+ Kz.) sin? g (57) 

Ky, = Ksg cos* + (Kpg- Kg,) sin@ cos $+ K,, sin® g (58) 


Ky Kp, Cos*'9 PK, - K,.") sing cos P+ Ke sin® g ie) 
Ky, = Ksp cos* @ + (Kss- K,y,) sini cos @+ Ky. sin? g (60) 
By, = Bpp cos* B + (Bag + Bep) sin cos P+ Bog sin? g (61) 
Byy= Bgg cos*  - (Bygt Bop) sing cos f+ By, sin? g (62) 


Byy = Bep sin? J + (B,,- Bg.) cos f sinf- Bys cos* (63) 


Computational Procedure 


The general case of unsteady motion of a circular cir- 
cumferentially grooved journal bearing has now been analyzed 
and the differential forces and differential spring and damp- 
ing coefficients are in forms which can be solved numerically. 
A detailed procedure, using a high speed computer, is given in 
Appendix B. However, the general computational procedure is as 
follows. For a given bearing geometry (L, R, C, andA/) and 
given operating conditions (WU, Z, n and é): 


Gir 


Reis | ays 


»~ ris © #506 gob sat) 


1. Determine direction of flow (+ 8) and locate up- 
stream and downstream cavitation boundaries. 

2. For each 9, calculate the film width. 

3. For each 8, determine the differential forces and 
differential spring and damping constants. 

4, Integrate numerically to obtain the radial and tan- 
gential forces and spring and damping coefficients. 

5. Calculate the bearing load, load angle, and horizontal 
and vertical spring and damping coefficients. 


Simplified Analysis 
It should be noted that equations 37 through 44, and, 


consequently equations 57 through 64 are independent of 6,8, 
Qe, and Ag due to the fact that (see Figure XI): 


Since dK,,,, 
independent of % and é in the above mentioned equations. 
If now frictional effects are neglected on the assump- 


tion that they are small compared to squeeze effects, then ~ 


and dBy~ are functions of w, then w must also be 


and UJ) can be used interchangeably throughout the entire anal- 
sis by virtue of the choice of reference axis. Then the equa- 
tions for the spring and damping coefficients and the equations 
for the steady state forces Fis and le can be non-dimension- 
alized and expressed in terms of just three independent vari- 
ables; eccentricity ratio, n, pressure ratio, PR, and dimen-~ 
Sionless squeeze, /A*. Let: 
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Then the pertinent dimensionless equations are: 


= {PR <F2))"Cio4%mR cos 6)" 


WwW 
“ly 5] Co a ee ae ea al ede lice 6 
L \/ 3n/A* sins a (05) 
Ww dw _ (#2 1 (1 + ncos6)°(PR - 1) (66) 
i a6 L° 2(cos @- cos 6*) 6n.A*(cos 6 - cos 6*) 
Bees pera (mys eds OW, (oe a) Gcs.6 de (67) 
ah sie L° (1 + ncos 0)* L ¥ 
cP... Ww sin? ) W 

= = i’ rt Oe, OS ee a Nae es ee ae ee Ae 2 Ls 
Sin =? SAE (=) as a ee) nee (68) 
CK pn C J® wi, sin® cosé 


(1 + neos 6 = 3ncos 6)dé6 


(69) 


eel 2 a al + ncos @) 


dK, Ae Wis sin? 8 al WW 7 
=| 7 (=) 9 ——________ + =(PR - 1) (= eld 
pin “| 2 (=) Tee Gees an ) (=) sin i SR Aig) 
dK. ¢ [ A® wy, Sin ®cos 6 1 | iw i 
a 6 ee nee ae lag 
yee) 
aKsy,c 3 WwW sin? ') 
== (Pp + Y¥ ore y 5 
mais Go wees eo) he) 
AB yy © cos* © 
ee MESS] gC Rear tinkchalon ame 5 
p Ik (1 + ncos 0)* oa (73) 
ABss ¢ sin® 6 
Bere Ae eee ett AG a 
p Lk /\ (1 + ncos 6)? (7m) 
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Cc sin8 cose 


ahs Ae 


E 
a; ae (1 + ncos 9)* 


ae (75) 


The general computational procedure is the same as before, 
except that all frictional effects are neglected. 
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Figures XIII, XIV, XV, XVI and XVII illustrate the shape 
and extent of the cavitation zone for various combinations of 
journal motion obtained by this investigation. Results obtained 
by Ishii and Richardson ( the latter in dashed lines) are also 
presented for comparison. The direction of the circumferential 
flow in Figure XVI is the reverse of the circumferential flow 
in Figure XV, thus Figure XVI is a mirror image of Figure XV. 

Tables 1 and 2 show a breakdown of the effects of circun- 
ferential squeeze, radial squeeze, feed pressure, complete film 
friction, and cavitation zone friction on the forces and spring 
and damping coefficients for the examples in Figures XV and XIV. 

Tabulated data for the simplified analysis which neglected 
frictional effcts is given in Appendix C for numerous combina-~ 
tions of dimensionless squeeze, pressure ratio, and eccentrici- 
ty ratio. The parameters were chosen so as to cover as wide a 
range of bearing applications as possible. 

| Table 3 iS a comparison between the data in Appendix C 
and data obtained by Holmes | 8] who analyzed the vibration of 

rigid shaft by extending Ockvirk's analysis and boundary con- 
ditions. The data in Appendix C was converted to Holmes nomen- 
clature as follows: 


peLE 
ra) QO 
een = Kyiy—— 
2Pi nm 2F 
ie Z f D, LR 
germ PORE 
Where: 
1. 
Ft hee (L-n*) +16n? IE 
PLR 4(1 -n?)* 
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which is a form of Ockvirk's load number. The factor of 
2is necessary because this investigation uses a circun- 
ferentially grooved bearing with a total length of 2L. 


Since Ockvirk's theory defines the attitude angle for a given 
eccentricity ratio, it was necessary to fix both @ and n in 
the comparison. 

The program in Appendix:B is written in a manner so that 
any or all of the values can be called without any major changes. 
Frictional effects can be neglected simply by deleting the in- 
Gicated cards, thus reducing the progr to three dimensionless 
variables as in Appendix C. It should be noted that for sim- 
plicity only one half of the bearing was considered in the 
analysis, but that the program and tabulated data are in terms 
of the entire bearing length 2L. 
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DISCUSSION OF RESULTS 


For the steady state example of Figure XIII, Ishii'ts results 
are very close to the results of this investigation, consider- 
ing that while Ishii used continuity of.flow through the cavity, 
he averaged the flow leaving the boundary l - 2 and then used 
this average to determine the boundary 2 - 3. Thus it should 
be expected that Ishii'ts point C should be outside the bound- 
ary of: this investigation since the average flow at Wo is less 
than’the straight line flow, while point D should be inside the 
boundary of this investigation since the avearage flow at Wp 
is greater’ than the straight line flow. However, h does not 
vary appreciably between point 1 and 2 due to the large slope 
of the upstream boundary, and thus the average flow through- 
out the cavity is very close to the exact flow. Also, some 
minor disagreement could be expected between machine and hand 
calculations. 

The deviation between this theory and Richardson's is 
quite pronounced along the downstream boundary, particularly 
in the steady state cases of Figure XIII and XIV; closer agree- 
ment is reached as the ratio of circumferential flow to axial 
flow is reduced as in Figure XV and XVI until there is perfect 
agreement in Figure XVII where only axial flow due to radial 
velocity é is present. None of the sample runs indicated any 
differences along the upstream boundary; this was to be ex- 
pected since both theories use the same method for calculating 
the upstream boundary. The values for load and attitude angle 
are surprisingly close in Figure XVII, considering that dif- 
ferent numerical integration techniques were employed for the 
two theories. Thus any deviations in the results can be attri- 
buted directly. to the differences between the downstream cavi- 
tation boundaries. This is best illustrated in the steady state 
examples of Figures XIII and XIV, where the most pronounced 
difference is in the fourth quadrantii(1l.5n1<98<2n)., 


46 


|= aA) ® 4 a q 
tal 4o° he! jade ‘ “ee fo) ite arr oe S14 


Eisaivatiiag (Winbiwicd cinttganwhd eis anco.6 O98 
Pee Sokol »¥ih few Tih eee Mh esses os4t4 Yreeky 
a2 alt inieresr tawetié Yo oient 40) 6a -aheaet @e 
Beas «te ti (ico EVE Oke Va weepss s- be equ Dos. 
mes 6s of C9{3 L4rts Gabe etude TLVe stens® al 
. 700 bit Pere. alee, Lies te eact .inseoty af a 
eeGs kaw sine PopRaged aaeerege ocd gnats. 
tae. yo belisee hapa a0 aetroads cted somdm 


@ 


Richardson's theory yields a greater film extent in this area 
and consequently a greater film force in the direction of the 
applied load, thus reducing thé net bearing load and increasing the 
load angle. -The same analysis can be applied to the other 
examples although it is somewhat more complicated in the dyna- 
mic cases. Richardson did not compute the spring and damping 
coefficients and so no direct comparison of these can be made, 
but by comparing the bearing loads and load angles for the 
dynamic cases;, it can be deduced that the same order of mag- 
nitude of differences in the spring and damping coefficients 
would exist by virtue of equations 45 and 46, 

The purpose of Tables 1 and 2 is to show the small effect 
that friction has on the film forces and spring and damping co- 
efficients, and therefore to justify neglecting friction in the 
preparation of the dimensionless tables in Appendix C. In gen+ 
eral, friction can be neglected in all cases except when 
(UW- 2%) and é are very small, as in half speed whirl. 

Table. 3 illustrates the: effect of the choice of cavitation 
boundary conditions on the spring and damping coefficients. 
Using Ockvirk's boundary conditions, Holmes assumes that the 
entire half of the bearing is cavitated , an assumption which 
approaches this theory only at the higher eccentricity ratios. 
Thus there is close -agreement at n = .6, while at n = .2 the 
two theories yield somewhat different values. At first, it 
would be expected that this theory should always give less 


values of Kr? K ib. wy rand B ng at low eccentricity ratios 


since these ey ee go to zero in the completely uncavi- 
tated situation(see Appendix C), but this is not so because 
the film and the film forces obtained by this theory are un- 
Symmetrical about 0=1.5 , and thus a clear comparison can- 
not be made. However, by comparing the results of this inves- 
tigation for a high and a low feed pressure, it can be seen 
that as the feed pressure increases causing the film extent to 


increase, the radial spring coefficient decreases. 
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Sternlicht( 13] used Sommerfeld's boundary conditions and 
digital methods to solve Reynolds! three dimensional equation 
and to obtain spring and damping coefficients for a full jour- 
nal bearing. However, no quantitative basis for comparison 
could be found, although it is significant to note that Stern- 
licht found the damping cross terms Bay and BL. negligible, 
while this theory and Holmes’ analysis found these coefficients 
to be of the same order of magnitude as Boy and Box 

The general validity of the results can best be ascertained 
by examining the initial assumptions that were made. Since 
low to moderately loaded bearings are considered in this investi- 
gation, bearing distortion and viscosity variations due to pres- 
sure can be neglected. Also the choice of a circumferential 
groove feeding arrangement reduces the variation of viscosity 
due to temperature since the oil flow for a line source is 
three to four times that for other feeding arrangements at any 
given load.14 . The accuracy of the short bearing approximation 
is illustrated in Appendix A where the error is negligible for 
L/D = +, and is acceptable for L/D = 4. 

The assumptions that undoubdtedly need considerably more 
investigation are those that apply to the flow of the fluid 
through the cavity,under dynamic conditions, particulserly 
where the striations join the downstream boundary. The steady 
state investigations(referenced in the Introduction) using 
transparent bearings have shown a finite and irregular distri- 
bution of striations, with the result that considerable mixing 
must occur at the downstream boundary, and thus in actuality 
the downstream boundary cannot be precisely defined mathema- 
tically. Also, because of the low pressure gradients immedi- 
ately surrounding the boundary, surface tension and wettability 
probably influence the shape of the boundary.: 

The other assumptions that were made such as laminar flow, 
Newtonian fluid, etc., are generally accepted and will not be 
discussed in this paper. 
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CONCLUSIONS AND RECOMMENDATIONS 


This investigation has analyzed and programmed the un- 
steady motion of a circular, circumferentially grooved journal 
bearing. The results should prove useful in designing the 
particular type of bearing considered, while the analytical 
and computational methods employed can be modified for appli- 
cation to other bearing configurations. The most important 
point that was brought out by this investigation was that. ne- 
glecting the oil film between 9 =n and 6 = 2n will not 
yield good values for the spring and damping coefficients in 
lightly loaded high speed bearings, where Ockvirk's theory 
that a given load establishes the attitude angle is invalid. 
Conversely, it would be expected that assuming a complete 
film in that region would also give poor results with the 
exception of hydrostatic bearings. 

Recommendations for future work in the subject area of 
this investigation are listed below. 


1. Devise and construct an experimental apparatus to 
measure dynamic forces and load angles, with sufficient 
accuracy to establish the validity of the boundary 
conditions of this investigation as compared to 
Ockvirk's boundary conditions. 


2. Extend Richardson's theory and his program to obtain 
the spring and damping coefficients. If the differ~- 
ences are relatively small, Richardson's analysis 
would be preferable due to the simplicity of his 
boundary conditions; in its general form it could 
more readily be modified for various particular anal- 
ysis such as half speed whirl, stability, vibrations, 
etc. 


3. For use with L/D>s, program Reynolds! three dimensional 
equation using the boundary conditions of contiuity. 
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APPENDIX A 
Accuracy of the Short Bearing Approximation 


Donaldson[?] obtained an exact solution to Reynolds' 
equation for a journal bearing using series solutions and ap- 
plying Sommerfeld's boundary conditions, but he assumed that 
the feed pressure was sufficiently high to suppress cavita- 
tion. ..Figure XVIII is a comparison of Donaldson's solution 
and the results of this investigation, using data for uncavi- 
tated bearings from Appendix C. Two length to diameter ratios 
are plotted as a function of eccentricity ratio versus the 
ratio of approximate to exact load capacity. 


Figure XVIII 


Approximate Load 
Exact Load 


fe) One 0.4 0.6 O58 a0 
Eccentricity Ratio 
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APPENDIX B 
Computational Procedure 


The analysis was programmed in Fortran II for use on the 
IBM 7094 digital computer. The program in its present form can 
be used either as a main program or as a subroutine by the in- 
sertion of the proper administrative and format statements. 

A sufficient number of comment statements is in tne program 

to permit someone who is familiar with the language to follow 
the logic, but some of the detailed computational methods need 
further explanation. 

It can be seen in Figure VIII that as the circumferential 
flow becomes much less than the axial flow that the downstream 
boundary conditions approach the upstream boundary conditions. 
and equation 22 applies throughout. If the flow is negative, 
Re@s, CU = 28) is negative, then the positions of the upstream 
and downstream boundaries are reversed. To determine which 
condition prevailed, the following expression was used: 


Velocity matte (VR) = circumferential flow oy Ane 
axial flow | Anc| +| é| 


This gave not only the direction of the circumferential flow, 
but also its relative magnitude. If VR was positive, all 
iterations were done in the positive 6 direction, and if VR 
was negative, then all iterations were done in the negative 
© direction. 

The beginning of the upstream boundary was obtained by 
iterating in 3 degree intervals, starting at 6 = m7 and calcu- 
lating w by equation 22, until the following conditions were 
Re’: 


veh Wy he) ge en Ls 
22. w(I) < w(I -1) 
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If these conditions were not satisfied when the program had 
iterated through 2n, then a non-cavitated solution was obtained 
by means of definite integrals (see Appendix D). 

If cavitation was present, then ©) was located by iterat- 
ing in one degree intervals starting at 6, , again calculating 
w by equation 22, until w(I) > w(I-1) i.e., the sign of the 
Slope of the cavitation boundary changed. This established 62, 
One degree intervals were used between 9, and 9, so that greater 
accuracy could be obtained in calculating h* later in the pro- 
gram. 

Then, starting again at ©, or 6, +3° whichever was even, 
the program was iterated in 6 degree intervals using equation 
22 to calculate the film width along the upstream boundary, and 
equation 27 to calculate the film width along the downstream 
boundary. Equation 27 was applied as follows: 


1. Assume that w(I) = w(I-1) 


2. Determine h* by iterating through the w's previously 
calculated between 90, and 6, , until w* > w(I). 
then interpolate between the stored values of ©* on 
either side of w(I) to obtain h*. 


dw 


3. Calculate the slope FG) 


h*(I-1) and w(I-1). 


by means of equation 27, using 


4. Then, w(I) ~ w(I-L) + oe Ae 


However, because equation 27 becomes invalid, both mathemati- 
Cally and analytically, as A--0, and because the quantity 
(h=h*) becomes very small and practically impossible to calcu- 
late digitally with any accuracy as VR — 0, it was necessary 
to apply certain restrictions along the downstream boundary. 
These were: 
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| >.25, equation 27 applied. 
|<.25, equation 22 applied. 
The number .25 was obtained by trial and error; as VR app- 
roached this value, the downstream boundaries by equations 22 
and 27 agreed very closely (Figures XV and XVI have VR values 
of .33 and -.33 respectively). At lower values of VR, this 
program gave inconsistant results for the reasons mentioned. 

For each © , the differential forces and coefficients 
were calculated by equations 28 and 28 and equations 37 through 
44, The integral of h*dz for use in calculating the friction- 
al forces in the cavitated zone was evaluated by Simpson's Rule 
using 3 ordinates 


e 
[xs ~ 2 (h* = h*)(w* - Wee + w) 
w 

Then the differential forces and coefficients were integrat~- 
ed from @=T1 to 9O= 3m, by means of the subroutine,which 
uses Simpson's Rule and 61 ordinates. For clarity, each compon- 
ent of each equation was calculated separately; thus DCSS1 is 
the first term of qKos 

‘The remainder of the program is self explanatory. It 
should be noted though that FR is the radial force directed 
inward; this was changed for conventions sake. 
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FORTRAN 


ANALYSIS OF A CIRCUMFERENTIALLY GROOVED JOURNAL BEARING 


INPUT IS TOTAL BEARING LENGTH, JOURNAL RADIUS,RADIAL CLEARANCE, 
ANGULAR SPEED;ORBITAL SPEEC,RADIAL VELOCITY,ECCENTRICITY RATIO, 
ABSOLUTE FEED PRESSURE,AND VISCOSITY (EL,RyCyOM,PHIDOT,DE,EN,PO,U) 
UNITS ARE LB,IN;,SEC,RADIANS 


CUTPUT IS RADIAL, TANGENTIAL, HORIZONTAL,AND VERTICAL SEARING FORCES 
SPRING COEFFFICIENTS,AND DAMPING COEFFICIENTS(FRyFT+FTOT,ANG ,CRR, 
CSS,CSR,CRSyGRRyBSS, URS, CXXyCVYsCY¥XsCXV ae BXXrBYYV a EX 


CIMENSTON X2(99),W2(99),X(61),W(61),DF11(61) ,DF12(61),DF13(61), 
10F14(61),0F22(61),0F23(61)»DF24(61),DCRR1(61),O0CRR2(61),DCSR1(61),% 
2CCSR2(61) ,AHSDZ(61),CCRR3(61),DCRS4(61)1D0CSS4(61),D0CSR3(61);5 
Sere othid y,Or25 (061) .W1 (124) 


COMBINE AND REDUCE INPUT 


GA=O0M-2.0*PHIDOT 
P=P0-15.0 
PI=3.1416 

FEE = 0.5*EL 

GG = (0.5%EL) #3 
©C=(U#R)/(2.0#C #3) 
CD=(UxOM#R¥#2)/C 
SS=DD/C 

TV=CC#OM 


BS=0. 
EW=0. 
PP=0. 
Be=0. 
CR=0. 
CS=0. 
CT=0. 
CU=0. 
CVv=0. 
CY=0. 
AHSDZ=0. 


MAGNITUDE AND DIRECTION OF CIRCUMFERENTIAL FLOW 


VR=GA/ (ABSF(GA) +ABSF(DE/(EN#C))) 
IF(VR) 2102110,110 


BOCATE X1 FOR NEGATIVE FLOW BY IFERATING IN 3 DEG INTERVALS IN 
NEGATIVE THETA DIRECTION, BEGINNING AT 5 PI 


CO 220 K=1,124 

PEL K= 124). 2115300.999 

A=K 

K=K 

X1=5.*P1+.05235—-(A#*.05235) . 
AH1LS=C#(1.0+EN*COSF(X1)) 
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AAA=COSF(X1) 
EBB=SINF(X1) 


EQUATION 20 


PH=.18 75#U#EL #%# 2%#(-2,.08#DE#AAA+GA# EN#C#BBB)/(Ce(LLO+EN#AAA) ) #%3 
TF(4SO#PH+P) 212,215,215 


EQUATION 22 


212 W1(K)=SQRTF(-((C#( LL O+EN#AAA) ) &#*#2eP)/(3.08%U% ((-2.0*#DE* 
LAAA)+(GA#EN#C#BBSB)))) 


iis 1) A999, 2204213 
Bae Pe OWL CK~=1)—-.499#EL) 2205214,214 
214 IF(W1(K)-W1(K-1)) 240,220,220 
215 wW1l(K)=0.5*#EL 
220 CONTINUE 


POCATE “X2°*FOR NEGATIVE FLOW BY ITERATING IN ONE DEG INTERVALS IN 
NEGATIVE THETA DIRECTION , BEGINNING AT Xl . 


240 


/ 241 
J)=X1-(A*®.01745) 

EQUATICN 22 
OW2(J)=SQRTF(-((C#(1.0+EN#COSF(X2(5))))#*#32#P)/(3.0%U#((-2.0#DE * 
BeEUSFUXZ(J)))+(GAXEN#CeSINFUX2(5)))))) 

IF(J-3) 250,250,245 

eae FF(W2(J)-W2(J5-3)) 250,251,251 
250 CONTINUE 
8251) WW2=W2(J-3) 

XX2=X2(J-3) 
JJ=J-12 
I=61-K/2 
eo.TO 153 


BOCATE .XLVFORIPGSITIVE FEOW BY ITERATING IN 3 DEG INTERVALS IN 
POSITIVE THETA DIRECTION,BEGINNING AT 3 PI 


meee £20 \K=1,124 
PRU EK= 124.21 by 300,999 

bal” A=K 
K=K 
KX1=3.0363+(A%.05235) +2.0#PI 
AH1S=C#(i.0+EN#COSF(X1)) 
AAA=COSF(X1) 
EBB=SINF(X1) 


EQUATION 20 


PH=.1875*UsEL #e 2% (-2.02DE#AAA+GA® EN*®C*#BBB)/(C#(LLO+EN#AAA) ) # #3 
PetsORP te) bo 2,059 15 é 
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EQUATION 22 


WL(K)=SQRTF(-((C#(1LSOFEN#AAA) ) ##3%P)/(3.0%U% ((-2.0*DE* 
LAAA)+(GA®EN#C#BBE)))) 


Hae) IOs P2051 13 
IFCWL(K-1i)-.499*EL) 12031142114 
IF (W1(K)—-W1(K-1)) 140,120,120 
WlL(K)=O.5#EL 

CONTINUE 


LOCATE X2 FOR POSITIVE FLOW BY ITERATING IN ONE DEG INTERVALS IN 
POSITIVE THETA DIRECTION , BEGINNING AT Xl 


BO, 150, J=1.99 

J=J 

A=J 
X2(5S)=X1+(A#.01 745) 


EQUATION 22 


OWw2(JV=SQRTF(~-((C#(1.0+EN#COSF(X2( J) 
) 


( Pte SEP VSS. ORUEC(—2e08DE) = 
ICOSF(X2(J))) +(GA®EN#*CH#SINF(X2(5)))))) 


) 
) 
Meh I =-3 Jo dS 09150145 
PeCW2( Jo —-W203-3)) 150,151,151 
CONTINUE 

hW2=W2(J-3) 

XX2=X2(J-3) 

JJ=J-12 

l=K/2 

Pete fIl=-K) to Ze) 559999 

I=[+1l 


INCREMENT IN 6 DEGREE INTERVALS AROUND THE JOURNAL BEGINNING AT 
mii, AL, EVEN) OR AT X1+3(IF Xi ODD) 


Boot yO N=1,61 

A=1 
X(1)=3.0#PI-.1047+(A#.1047) 
AAA=COSF(X(T)) 
EBB=SINF(X(I)) 


IF RELATIVE MAGNITUDE OF CIRCUMFERENTIAL FLOW IS SMALL, 
EQUATION 22 APPLYS THROUGHCUT 


Prec iciotst VRy— ido)" bed 594. 159 
IF(VR) 2604,1604,1604% 


IF (XX2-X 


GU) a b54 Wl 540166 
Binne—X{1)) 166,154,154 
EQUATION 20 


PH=.1875#U#EL ## 2%(-2.0*DE*AAA+GA*® EN#C*#BBB)/(C#¥(L.OTEN#AAA) ) % #3 
LEC 4.0%PH+P)° 165,1675167 
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UPSTREAM W (EQUATION 22) 


165 WC IT) =SQRTF(-((C#(1.0+EN#AAA) ) #23 2P)/(3.08Ue ((-2.0#DE#AAA) +(GAHEN 
1 *C*BBB)))) . 
AHS=C#(1.0+EN#AAA) 

XS=X(1) 

GO TO 168 


COWNSTREAM W IS COMPUTED FROM STATEMENT 166 THROUGH 1667 


eee ke (VR) 11515315 
ELS R Perl =6 lL ) 12,13,999 
12 WOT )=WC(I+1) 
GO TO 20 
13 w(T)=WC1L) 
so TO 20 
15 Pa STs bg 7, 
16 w(I)=W(61) 
eo TG 20 
Md wW(T)=W(I-1) 


He FOR EQUATION 27 


(20 CO 1664 J=1,JJ 
fers) "S67 16 5,999 
36 PF(W(T)-W2(5)) 1664,1661,1661 
mooie IFC J—1)999,1662,1663 
1662 XS=X2(J) 
Ge 10 1665 
Geos xs—X2(J—1 CW (1 )-W2 05-1 de CX 205 X20 5-1) ) J CW 209) -W2 05-1) ) 
moet: L665 
1664 CONTINUE 
1665 AHS=C#(1.0+EN®COSF(XS)) 
AH=C*(1.0+EN*AAA) 


EQUATION 27 


CwW=(0.5*#w(1)#(-2.0#DE#AAA+GA#EN#C *#BBB)+( (AH®*3%P)/(6.0#UFW(I)))) J 
1(GA#(AH-AHS) ) 
1667 wW(IT)=W(1)+ABSF(OW)#.1047 


59 PeOW(l }-EL*#O.5) 1685-167 ,167 
meot = W(l)=0.5*EL 

©1168 AA=(1.0+EN#AAA) &%3 

BB=W(I) #23 
CQ=(1-0+EN#AAA) #82 

RR=QQ# #2 

YW=(1.0+EN#AAA) ##2 


CIFFERENTIAL FUNCTIONS OF FILM WIDTH AND TRIGONOMETRIC FUNCTIONS 
IN CAVITATED ZONE 


memuve ALL CARDS DOWN 10 1683 IF FRICTIONAL EFFECTS ARE «TO BE 
NEGLECTED. 


1682 AHSDZ=0. 
| IF(.49*#EL-W(1)) 1683,1683,1671 
1671 AH1L2S=AH1S+(AHS-AHLS)/2.0 
a COD=((AH12S/C)-1-0)/EN 
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X12=ACOSF(O00) 
EEE=SINF (X12) 


EQUATION 20 


PH=.1875%U#EL #%# 2% (-2.02DE#COD+GA# EN#CHXEEE)/(C#(1.0+EN*DDD) ) #¥3 
IF(4.0®PH+P) 1673,1672,1672 


Poe rE E=—E GE 
EQUATION 22 


1673 wWi2=SORTF(—((C#(1.0+EN#® 00D) ) ##3%P)/(3.08U%((-2.0*DE*DDD)+(GA*EN 
ie GY EEE) ))) ) 


EQUATION 22 


WS=SORTF(-((C#(1L.-O0+EN#COSF(XS) ) )##3*#P)/(3.0%Ue( (-2.0*DE* 
BCOSF (CXS) ) +(GA#EN&C#SINF(XS))))) 


SSS=ABSF ( (AHIS—-AHS)/2.0) 
AHSOZ =(SSS/3 .0)*#(W1(K)+4.0*W12+WS) 


1683 CCRR3(1)=AHSDZ*®*BBExAAA/AA 
CCRS4(1) =AHSDZ#AAA/QC 
CCSS4(1)=AHSDZ*BBB/QQ 
CCSR3(1)=AHSDZ*#AAAR#2/AA 
CF15(1)=AHSDZ*B8B6/QQ 
CF25(1)=AHSDZ#AAA/QO 


CIFFERENTIAL FUNCTIONS OF FILM WIDTH AND TRIGONOMETRIC FUNCTIONS 
ae COMPLETE FILM ZONE 


CF1L(I)=(BBeAAA%#2)/AA 
CFL2(1)=(BB*AAA*BBB)/AA 
CFI3(I =W(I)#AAA 
| CFL4(1)=W(1) *BBE/(1L.O+EN#AAA) 
CF22(1)=(8B*BBB*#2)/AA 
CF23(1)=wll)=88B 
CF24(1)=W(1) *AAA/(1.0+EN#AAA) 


| _ EPCRR1(1)=B8B* BBB #AAA¥#2/RR 
CCRR2(1)=W(1) #AAA*BBB/QQ 
CCSR1( 1) =B8B*8BE#*2#AAA/RR 
CCSR2(1)=W(T) #AAA##2/QQ 
IF(.49#EL-W(1)) 1685,1685,420 

420 N=174+6«1 
1F(M-360) 422,422,421 

421 N=M-360 

422 ZZ=2.0#W(1)/EL 


Miser RaNT SSTAVEMENT HERE “PF OTRETA AND FILM WIDTH(PER CENT) IN 
BAVITATED, REGION TS. DESIREC. 


1685 IF(VR) 1686,1688,1688 
1686 1=I-1 

IF(I) 999,1687,i70 
1687 1=61 

Bit Ou 200 
1688 I=I+1 

PEL -62)— 110, 1689, 999 
1689 1=1 
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CONTINUE 


INTEGRATE FUNCTICNS OF FILM WIDTH AND DIFFERENTIAL TRIGONOMETRIC 


FUNCTIONS,USING 


PA =SIMP(DF11) 
BC y= SI MP COFL2} 
BD = SIMP(DF13) 
BG S= STMP{ORZ2) 
Bio JS (MP DEZS) 


ER=SIMP(DCRR1) 
EV=SIMP(DCSR1) 


SeeE LE NEAT. of 


BW=SIMP(OCSR2) 
BS=SIMP(CCRR2) 
EE = SIMP(DF14) 
EP = SIMP(DF24) 
R=SIMP(DCRR3) 
CS=SIMP(DCRS4) 
ET=SIMP(DCSS4) 
CU=SIMP(DCSR3) 
EV=SIMP(DF1i5) 

CY=SIMP(DF25) 


COMPONENTS OF 


Fll=-2.0+#*CC#DE« 
F12=CC*#GA#EN*C# 
F13=0.5#P#R#ED 
Fl14=-DD#BE 
Fl15=-SS#DV 

hele —2 .O*CCeDEX 
F22=CC#GA#EN*C% 
F23=0.5*P#R*BH 
F24=DD#8P 
F25=SS«DY 
Sanl=11#(B8C-3.0 
CRR2=SS#BS 
CRR3=2.0*SS#DR/ 
BRR=-2.0#CC#BA 
CRS1L=-TT*BG 
GAS2Z2=-P*#R*BH/2. 
Gaoo=—-SsS*BP/EN 
CRS4=-SS#DS/C#E 
BPao=—2 .0+CC#BE 
Sosl=Ti*eBc 
CSS2=P*R*BD/2.0 
Bosoo=—SS#SE/EN 
So54=-SS*DI/C#xE 
BSS=-2.0*CC«*8G 
CSR1=TT# (BG-3.0 
CSR2=-SS#BYW 
CSR3=-2.0#SS*DU 
BSR=-2.0«*CC#BC 


BUNT O. 350 


SIMPSONS RULE 


CARO Sm. bs-FRLECELGNAL ERFECTS ARE 10. BE 


EQUATIONS 28,29;AND 37-44 


BA 
BC 


BC 
BG 


*EN#BR ) 


C 


OxCxEN 
N 

*C#EN 

N 
e#EN*BY) 


/C 
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NEGLECTED 


32,0 


B39 


FORCES AND SPRING ANO CAMPING CONSTANTS IF NO CAVITATION 


PRINT 912 

AL=-PI#EN/(1.0-EN##2)##2.5 
AG=(3.14/EN)#(1.0-(1.0/(1.0-ENe #2) ##0.5)) 
AC = 1.57/(1.0-EN##2) 481.5 

Fl1l=0. 

Fl12=-0.0 

F13 2-00 

Fl4=0. 

F15=0. 

F21=0. 
F22=CC#GA*GGKEN#2 .0#AC&C 
F23=0.0 

F24=0. 

F25=0. 

CRR1=0. 

CRR2=0. 

CRR3=0.0 

BRR=0.0 
CRS1=-GG*TT#2.0#AC 
CRS2=0.0 
CRS3=-(SS#EE/EN )#2.0%AG 
CRS4=0.0 

BRS=0.0 

C551=0..0 

6ss2=0..0 

CSS3=0.0 

CSS4=0.0 
BSS=-GG#4.0*CC#AC 
CSRL=GGeTT#(AC-3. 0% EN#AT ) 
CSR2=0.0 

CSR3=0.0 

ESR=0.0 


baa OX (Eero tr Late iS) J 7. Set beerR) 
ri=7.O0# (VR 2itpZe t+F2atF 244725) 7 UT «SHEL ER) 
PrOtT=SQORTE CE Tes 2rrR ee 2 } 

PHI=ATANF(FT/FR) 

ANG=57.3#PHI 

CRR=—2 OH (CRRIFECRRZ+CRRE )#C7( 7.5 EL eR) 
Po=—2 «O0* (CRSIFERSZECRS3S4CRS4 V*¥C/J(7.5*EL eR) 
Pos —2 00 e( CSSlresszZ+CSS37CSS4  7*C/ (7 s5*EL *R) 
mon—2.O* (CSRIFCSRZFESRS )*#C/(7.5* ELLER) 
BPBRR=-2.O#BRREOM*#C/( 7.58 EL#R) 
PSs=—2.0*0M#C/(7.5#EL*R) 
BRo=2.O*BRS*OM#C/ C76 5SeEL BR) 

AAA=COSF (PHI) #*2 

EBB=SINF( PHI) ##2 

CCC=COSF(PHI)*SINF(PHI) 
CXX=CRR*AAA+(CSR-CRS) *CCC+CSS*BBB 

ET Y-Grk*oBb+t(CRS-CSKISCCEFCSS AAA 
CAY=CRS*AAA+( CRR-CSS ) #€CC+CSR*BBB 
CYX=CRS#BBB+(CSS-CRR) «CCC+CSR#AAA 
EXX=BRR*eAAA+2.0%BRS#CCC+BSS*#BBB 
BYY=BRR#BBB—-2.0#BRS#CCC+BSS#AAA 
EXY=BRR#CCC—BRS «AAA-BSS*#CCC+BRS#*BBB 
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FORTRAN : 
SIMPSONS RULE FOR 61 ORDINATES, WITH 6 DEGREE INTERVALS 


FUNCTION SIMP(Y) 

CIMENSION Y(61) 

COD=0. , 

a EVEN=0. 

ie £0*800" 1=2,60,2 

 EVEN=EVEN+Y(1) 

CO 810 1=3;59,2 

y~ CDD=ODD+Y(T) , 

— SIMP=(.1047/3.0)#(Y(1)+4.0#EVEN+2.0#00D+Y(61)) 
RETURN 

~* END 

ye CATA; . —< 
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APPENDIX C 
Tabulated Values of Bearing Load, Attitude Angle, 
And Spring and Damping Constants 


Tables 4 through 45 contain values of dimensionless load, 
load angle, and dimensionless spring and damping constants as a 
function of eccentricity ratio, pressure ratio, and dimension- 
less squeeze. The eccentricity is in increments of 0.1 from 
0.1 to 0.7, the latter being chosen as the upper limit in 
order to remain within the limits of the short bearing approx- 
imation, the assumption of constant viscosity, and the assump- 
tion of no journal or bearing distortion. The pressure ratios 
are 1.5, 2.0, 3.0, 4.0, 5.0, and 6.0, which adequately cover 
the range of most bearing applications. The choice for varia- 
tion of. squeeze was made as follows; from Shaw and Macks 


; 


-003 


I 
IN 


small bearings: 


~ 
; 


~0005 


a 


fe fee 
v 


large bearings: 


Now let -- = Vary from = to 1. Then R varies from 2L to $L 


and: 


z 
ois 3 f ab iz 5 
/ i] =| = 278) x: 107 
C | 1 2(.003) 

L Le L 3) 
Peele ois 72 

om 2 7 

awe =! | ad 1.6 axe 10 
Eas | 0005 | 


Now assume a viscosity range from 10° to LOL. lb-sec/in”? and 
@ speed range from 100 to 2000 rad/sec. Then: 


/A* = = (hy 


Cc 
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ms *cbocl 


= % 
€06. j94 1". 


Ax = 4405°) (2.00) (2780207) Lage 


MIN 15 


(107° )(2000)(1.6)(10’ ) 
SN EN epee € ey SR eae 


However, no bearing applications at the upper limit could be 
conceived and so it was decided to vary /A\* from .1 to 8000 
in 30 steps that would permit easy interpolation. It should 
be noted that generally the first few entries in each table 
are uncavitated bearings, which can be recognized by virtues of 
the fact that = 90° 

The tabulations are in terms of the entire bearing length 
(2L) and are converted to dimensional values as follows: 


F=F p, LR 
c 
Km = a amas 
) 
Be = BNM 
Pa 


For end fed bearings of total length L, divide the above 
values by 2. 
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APPENDIX D 
Table of Integrals 
The following definite integrals are used in evaluating 


equations 28 and 29 and equations 37 through 44 for the non- 
cavitated case, 
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